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The ID and 2D Bose-condensation of trapped atoms in a gravitational field 
are considered. The deformation of the finite parabolic potential in this field is 
modeling via the combination of two rectangular ID and 2D traps, for which 
the cut-off and the re-definition of spectrum are taken into account. A Bose- 
condensation T c shift by the gravity is calculated. A sign and a magnitude of it 
in a deformed trap depends on the order of including the gravitation field. The 
special choice of this order may describe three consistent Bose-condensations 
with different temperatures. These transitions may be associated with a trans- 
portation of a trap on the cycle (I) Earth-(II) Space-(III) Earth. 

PACS: 05.70.Jk, 67.40.-w 



The news on the critical temperature T c of the Bose-Einstein condensation (BEC) 
entered the last decade. It is known, that the critical T c of the superfluid 4 He (the oldest 
- though inderect - example of BEC) is the same for any tube with a liquid helium. The 
decreasing of 4 He T c in porous glasses was found experimentaly [1] and was explained 
theoreticaly [2,3], and this shift of a critical temperature depends on the parameters of a 
porous glasses [4] . Meanwhile the calculation of BEC T c in an infinite volume for bosons 
with a weak interaction between them showes the increasing of T c [5,6]. 

As to the ideal gas, it was found [7], that in an external field, restricting an infinite 
volume with this gas, T c depends on the form of a field. The experimental discovery of 
alkali atoms BEC in magnetic traps [8,9] gives the opportunity to discuss the same prop- 
erties in finite systems. A finite volume assumes mesoscopic boundary conditions, which 
require a new formulation of the thermodynamical limit in theory [10]. The decreasing of 
an ideal gas BEC T c was shown [11] as a finite size correction N ~ oo — ► N < oo for the 
"trap version" thermodynamical limit hu/T — > 0, N — > oo calculations. 

Here we intend to look for a change of BEC critical temperature T c via the influence 
of a gravitational field on T c of a BEC gas in traps. A problem of the gravity influence 
on the critical temperature of 4 He arises in the connection of Space experiments on this 
subject, and the same experiments on BEC in trapped atomic gases are expected in the 
nearest future [12]. So it is interesting to predict theoretically the contribution of gravity 
field to the critical BEC temperature in traps. We deal with the ID and 2D theory of 
spin-less non-interacting atoms in a finite trap. The basic parameters are taken from the 
BEC of alkali atoms experiments [8,9]. We start with the isotropic parabolic trap with a 
potential barrier Uq, frequency u and a height h and introduce two rectangular traps in 
order to describe the main part of gravity effect via the deformation of a trap potential. 
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The potentials of a parabolic trap are different for the cases of a zero- and non-zero 
gravity field in direction z. The potential inside of a trap may be written as 



U(z) =U g + (mu 2 /2) (z + A) 2 , -h/2<z<h/2, A = g/u 2 , (1) 

where m is the mass of an atom, g is the free fall acceleration constant. The dashed line 
is an initial parabolic potetial, a solid line is a parabolic potential, deformed by a gravity 
field. The gravity field shifts the minimum of potential from its center z = by the 
value A and down by the value U g = —mg 2 /2uo 2 . The difference between new (shifted) 
potential barriers U + (on the right) and LL (on the left) 

2 / h \ 2 1 

u±= ™_l A ±_\ UQ _ U =U+ _ Uo = _ mgh 

is equal to mgh. It follows from these formulas that the condition of macroscopic stability 
of the trap in gravity field is g < u 2 h/2. It is fulfilled for the experiments [8,9] and the 
last inequality will be satisfied for the frequencies more than 100 Hz. For all that the 
quantum dynamics of atoms undergoes nonperturbate gravitational disturbance because 
of the big shift of the trap potential (U g ^> ftiv) in comparison with its frequency. 
One-particle boson eigenfunctions u n and their energies e n > for the 2D-case 

{2m ~) 
V 2 + [e n - U(r)]ju n = 0, r = {x,z}, u n (r) = u nx (x)u nz (z), (2) 

6 ™ = 2m / ^ Un ^ Undr + / U(r)u n u n dr, n = {n x ,n z } 

are the starting points of the quantum analisys of any finite trap U(r), a ground level of 
it is £o 7^ 0. In order to introduce the Gibbs statistics we use the path-integral device, 
applied earlier [13,14] for the BEC of interacting bosons. 

The system of iV bosons (2) in a volume V may be represented by the Hamilton 
function H 

v h 2 



H -> if = dr 



— (Vip*(r, t)Vip(r, t)) + i>*U(r)i> 



Here ip*, ip are the path-integral trajectories with periodical boundary conditions on [0, f3}. 
Following [15,16], we separate the trajectories 6q, bo as the "slow" variables relative to the 
"fast" trajectories 6* , b n , n ^ for T < T c 



r(r,t) = -L I J2 b* n (t)u n (r) + b* (t)u (r) ] , |6 | = N » 1, 



so that the function N (t) is the path-integral image of the " Bose-condensate" particles 
in Bogoliubov theory [15]. The "broken gauge symmetry" is known as a property of the 
bose-gas with the separated condensate fraction. Still the quasiclassical integral of motion 
iV = iVo + M for the system (2) may be expressed by the equation [13] 

^ = {H, N } + i[H, JVx] = 0, N 1 = J2Kb n - (3) 
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Here Aq is a number of non-condensate bosons, {, } is the Poisson bracket in 6 0) &o classical 
amplitudes and [, ] is the quantum commutator via operators &+, b n >, n, n' ^ 0. The action 
5" with the kinetic term K for iV non-interacting bosons in any volume V looks like 



dt. 



An effective action S e f(p,p) for a condensate density p = N /V and a chemical potential 
p is defined via the partition function with constraint 

Q = J e s <tM dp dp = Sp [exp(-(3H)5 RtRl+p ] , R = N/V, R 1 = NjV. 

The equation (3) is the background for the definition of p, and by the use of Fourier 
decomposition for 5r^r 1+p with p = iy/j3 we represent Q 

Q= I Dbl Db [] / Db* n Db n r dy exp[ty(R l + p - R) + S/h] 

as the integral over y, the functional integral over the "slow" bo,b^ an d "fast" functional 
integral trajectories b n ,b* n , with the periodical boundary conditions in [0,(3]. The calcula- 
tion of Q over "fast" variables gives the effective action S e f of the condensate 



Q = J DblDb J dpexpS e f, 



S, 



ef 



-b. 



db 
1 dt 



p(p - £ ) - i^R 



dt - In (1 - exp[-/3(e n - p)]) 



It leads to the variational equations 5S e f(bo,b1j, p) = in the form 



dS, 



ef 



(3dp 

f(tn x ,n z ,(3,p) = f 



p + -E- R = 0, E 



exp[P(e n - p)] - 1' 



n = {n x ,n z }, 



(4) 



dS, 



ef 



Pdb* 



dbo 
"dt 



bo(e - p) = -Lb = 0, p — e - kT\n ( 1 + 



iVn 



(5) 



Equation (4) means the balance of condensate and non-condensate particles. The solution 
of equation (5) follows via the calculation [13] 

(^-^[exp^eo-^-l] 

of the operator L^ 1 with the periodical boundary conditions. The equations (4,5) deter- 
mine an BEC solution for the model (2,3) below the critical temperature for the case of 
translation-noninvariant systems. So, the critical temperature T c for a finite-size trap is 
determined by the equations S(T C ) = N, p = 0, while at the temperarure T ~ the 
equations p ~ R, R\ ~ 0, p ~ e are present. Note, that if we take 6q, b as constants, we 
get only T = limit of solutions for p, it is a property of the ideal Bose-gas. The formulas 
(3-5) are valid for any trap U(r). We apply them for a parabolic trap in formulas (6-8). 
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The best way to catch a mesoscopic effect is to look for the shift of a critical tempera- 
ture T c of BEC in a trap. It is obvious, that a finite size h for a parabolic trap manifests 
itself in the following mesoscopic properties: 

(i) a finite number of levels n max = Uq/Tiuj of the initial symmetrical parabolic trap, 

(ii) a diminishing of the left-side potential barrier U — > (£/_ — U g ) < U by the gravi- 
tational field. The re-definition of one-particle matrix elements e n must be done in both 
these cases. It means that we must (i) cut-off an upper limit n < n max for barrier Uq 
and calculate energies e°, as well as (ii) cut-of an upper limit in (2) for barrier U- and 
calculate energies e*. So the further calculation will be done with the formula 

Umax ( p _ Q 

N = N + N 1 , N 1 ='£f n , e n = \ f 9 ', n = {n x ,n y }. (6) 

of the full number iV of atoms in a trapped ideal bose gas. The level of n max defines the 
upper bound states of atoms. The numerical parameters of the experiments [8,9] for a 
parabolic trap in an isotropic version may be reduced approximately to the values 

U ~ 10- 9 eV, h = 2mm, U + -U^= mgh, uj ~ l()- 13 eV. (7) 

We suppose the critical 3D peak density [8,9] R ~ 10 14 cm~ 3 to be independent on time 
as the evaporation cooling is neglected here. The nD scaling for the number of atoms in 
a parabolic trap with a volume V, over wich the system is confined, looks like [10] 

R\nD = V , V~l — \ , 71=1,2,3. (8) 

In order to symplify a problem we would like to represent the parabolic trap by the 
rectangular trap. It means, that we describe the deformation of the parabolic potential in 
the gravitation field by transforming the symmetrical rectangular trap with potential Uq to 
the asymmetrical one with barriers U + — U- = mgh. Our purpose is to estimate the shift of 
a critical temperature in the asymmetrical rectangular trap relative the symmetrical one 
as a measure of the gravity field influence on BEC, tunneling effects are not considered. 
The values of R in a ID and 2D parabolic traps are 10 5 cm _1 and 10 10 cm~ 2 in the view 
of (8). We keep these R the same for a rectangular traps. So a number of atoms in a ID 
rectangular trap is Ni\id ~ 10 5 cm _1 x 0.2cm ~ 2 • 10 4 , and in a 2D rectangular trap is 
N\\2D ~ 10 10 cm~ 2 x (0.2cm) 2 ~ 4 • 10 8 . The equation for a spectrum e n — > Ek for the 
symmetrical ID and 2D traps looks like [17] 



/I1N 2kh 
tg(kh) = — 

Pk 



k 2 h 

{Pkf 



h 2 (k 2 + k 2 ) i 

E k = z \ p k = ^j2m(U -E k ), k = {k x ,k z } 



2m 

and for the asymetrical (deformed) trap for k = k z looks like 

k 2 Ti 2 




tg(kh)=k%[ — + — )[^ r ^-l) , pt = yj2m(U ± -E k ). 



The distances between the energy levels E k of an initial rectangular trap with the given 
parameters (7) are ~ lO~ l9 eV at a bottom and ~ 10~ 13 eV at a top barrier Uq. It has 
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~ 10 5 levels (contrary ~ 10 4 in experiments [8,9]), so the average energy distance between 
levels ~ 7 ■ 10~ 1A eV is the same as the frequency huo ~ 6 • lO^^eV of the paprabolic trap. 

In order to get more systematic picture we start with ID case with the R ~ 10 5 cm _1 
density. The critical temperature T®\ie> = 2.41 • lO" 10 -^ for an initial rectangular trap Uq 
is found from the equation (6) with the use of calculated spectrum E k 

No = 0, N°\ 1D = N\id = E f( T c\iD,El) = 2 ■ 10 4 , k < k max . 

A suitable accuracy of calculations may be received for k max ~ 1000 in the last sum in k z 
variable. The deformed rectangular trap U- ~ 10~ 8 t/o is determined via the parameters 
h and U + — U- = mgh in (7). It contains only ~ 10 levels, which are occupied by 

10 

N£\id = E f(Tc\iD,El) = 1.8 • 10 4 < N\ 1D 

atoms. These atoms undergoes a new BEC transition in the system of re-defined levels 
with a critical temperature in a deformed trap Tf- \\d = 2.27 ■ 10~ 10 , that is found from 
the equation 

10 

N = 0, N?\ 1D = E f(T?\iD, E£) = 1.8 • 10 4 , k < k max . 

So the shift of the critical temperature is A(T c )\id = T^Ild — T®\id = —0.14- 10~ 10 . The 
rest part 0.2 • 10 4 of atoms becomes non-trapped (continuous spectrum). 

As to the 2D case, the sum over the variables k x , k z is taken within the constraint 
condition [] E k < E(n\ max ), k = {k x ,k z }. A temperature T°| 2 _d = 5.05 ■ 10^ 7 i^ for an 
initial rectangular trap Uq is found from the equation (6) 

N = 0, N X \ 2D = N\ 2D = E f(T°\ 2D ,El kz ) =4-10 8 , k < k max . 

k x ,k z ^0 

The latter sum includes all the items of the partition function within the possible degen- 
eracy of levels. A suitable accuracy of calculations may be received for k max ~ 400 in the 
last sum in k x , k z variables. The deformed rectangular trap LL ~ 10 -8 t7o contains ~ 10 
levels, which are occupied by 

10 

N?\2D= E f(T?\ 2 D,E° kxk J = 0.96- 10 s <N\ 2D 

k x ,k z ^0 

atoms. Just like in ID case, we look for a transition temperature in a deformed trap using 
an equation 

10 

N = 0, N±\ 2D = E f{Tf\2D,El kz ). 

k x ,k z ^0 

The result of calculation is T^\ 2 £> = 5.03 • 10~ 7 K, so the shift of a critical temperature is 

A(T C )| 2D = Iflu, - T% D = -0.02 • IQ- 7 K. (9) 
lr This condition is trivial for an initial trap, but rather effective for a deformed trap. 
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Now we describe the T c shift in a system with the same number of atoms, 4 • 10 s for 
example (for 2D-case in our theory). The critical temperature for the symmetrical trap 
was found above as T c °| 2 _d = 5.05 ■ 10~ 7 K. The critical temperature for the asymmetrical 
trap (T±\ 2D )* = 2.02 • lO'^K is found from the equation 

10 

iVo = 0, ^{20 = ^20= £ f[(T?\2 D ) ,E± kz ] = 4-10 8 , 

A(T C ) = (T±\ 2D y - T C °| 2D = 1.52 • 10~ 6 K. (10) 

Note, that (T^\ 2 d)* in (10) does not coinside with the T^\ 2 d in (9), the reason is that 
the processes (9) and (10) are different. This important fact will be discussed in the end 
of the paper. 

The above calculations concern a rectangular trap. As to an initial parabolic trap, it 
has the same barriers U± with the same shift U+ — U- = mgh, but differs in spectrum. It is 
reasonable to suppose, that the same scale of a BEC critical temperature shift will be valid 
for the parabolic trap in the gravitation field, so the qualitative correspondence between 
the critical temperatures T c ° — > T\ g=0 , Tf — > T\ q ^ for a rectangular and parabolic traps 
is present. Meanwhile the tunneling effect is not reproduced by a rectangular trap. As 
it is seen from the formula (1), the macroscopic stability of a parabolic trap in a gravity 
field is expressed by the unequality uo 2 > 2g/h. It is clear, that this condition is trivial in 
the absence of gravity field or for the trap on the whole axe — 00 < z < 00. The finite size 
trap effect in the case of a deformed rectangular trap n max ~ 10, attributed to the most 
noticable gravity shift of T c , may be reproducrd for the parabolic trap via the condition 

(uo 2 h-2g) -> +0. 

This condition means, that the largest shift of T c may be noticed just before the destruc- 
tion of the parabolic trap. 

Let us turn in a conclusion to the non-ideal gas with an interaction G(r, r') beteen 
atoms. In the ID case the effective action in (4,5) for interacting bosons is generalized in 
the Hartree-Fock-Bogoliubov approximation as [13] 



^ nn'^O I 



S ef (p, /i) = -Pp 2 l0 V + (3p{p - R)V - Ppe + 

n 



P(W nn , - p5 nn ,) - Aln sh^- + PA 



F 2 



N 2 m= /-l 0\ M - ( W nn' - I^S nn ' 2^ nnl p 



V ' \ 1 y ' ' y -27n„'P -W nn > + p8 nn > y' 
W nn , = w nn /V + 2p 7nn ,, E 2 nn , = (M\ u ) 2 - (M| 12 ) 2 , <\? n = (b* , b )(p l0n ). 
Matrix elements of G(z, z') in ^-approximation G(z, z') — > G ■ S(z — z') 

G f , a G f , o G 



7o 



= — J dzuQ, 7on = — J dzu 3 u n , 7„ n / = — J dzu 2 u n u n > 



are detrmined for the functions of harmonic oscillator in z direction with the potential 
barrier (C/_ — U g ) < U at the left endpoint of the trap. Only the cut-off spectrum 
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contribution (without the re-definition of e n ) is taken into account. Equations SS e f(p, p) = 
for a trap with the constant density of atoms R = N/V may be evaluated approximately 
by dividing the sum over full number of trap levels into two terms with n<no and n ^ no. 
Parameter n ~ 100 is defined under the condition 

n hu ~ p7„ ono ~ GN , N ~ N, (11) 

following from the unequality 



10- n eV -hu<^ GNo ~ 10 _11 eV; G = 



m 



where the energy of interaction between atoms is taken from the estimation [8] for the 
scattering length a = 4.9nm with density of atoms R ~ 10 10 cm~ 3 . Thus, a density of 
condensate is written in the form p = p < + p > , where p < = p|„<„ , p > = p| n >n - We can 
estimate the ratio Y of Bose-condensate densities in a trap with p\ g and without gravity 
p|o 

Y _P\g _ (P < + P > )\g ^ P K \g (-, P > \o P > \ 9 \ <, >, 
P\o (P < +P > )\o P<|o V P lo P < l 9 / 

In the case of a "strong" gravity field (LL — U g ) C/o the inequality p < \ g > p < |o follows 
after the evaluations. In the case of a "weak" gravity field ([/_ — C/ fl ) ~ C/o the equations 
p < | 9 = p < | and p > \ g > p > \o are valid. The inequalities 

y -Plk(i-plk\>i y . - [\ p > \°~p > \A i 
s * rons p<| ^ p<| ; ' w ~ v p<| ; 

are proved. Taking into account, that the large Bose-condensate density corresponds 
the large critical temperature, we note, that the gravitional field increases the critical 
temperature T c of Bose-condensation for non-ideal gas in a trap. The parameter no in 
(11) corresponds the parameter n max in (6) as the border of BEC occupated states. The 
tunneling effects are also neglected here. 

We would like to look for a correspondence between formulas (9,10) and the probable 
motion of a trap between the Earth and Space. We measure the (Tj^d)* = 2.02 • 10~ 6 fT 
temperature for the N\2D atoms on the Earth in a deformed trap in the begining. Then the 
trap goes to Space, and the same number of atoms N\2D undergoes the phase transition (in 
a non-deformed trap) at the temperature T®\ 2 d = 5.05-10 _7 i^ with a shift (10). If the trap 
goes back to the Earth, some part of atoms becomes non-trapped due to its deformation. 
The rest part N^\ 2 d of them undergoes a phase transition at the temperature T^\ 2 d = 
5.03 • 10 -7 fr, and a temperature shift is given by (9). 

So we have three transition temperatures (with the same trap!) in different cases I, 
II, III 

(Tf\ 2D y = 2.02 • 10~ 6 K -> T C °| 2D = 5.05 • 10' 7 K -> T?\ 2D = 5.03 • lO" 7 ^, (12) 

start on the Earth (I) — > in Space(II) — > finish on the Earth (III). 

This is a picture for an ideal Bose-gas. The sign of T c shift for the non-ideal gas on 
the stage I— > II was performed in the last formulas of the paper, but the re-definition 
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of levels was not taken into account. It seems, that the 3D calculation will increase the 
magnitudes of the critical temperatures in (12). 

The mesoscopic nature of the problem considered above is seen in the strong depen- 
dence of the calculated T c shifts on the initial parameters, and the way of their association 
with experimental prototypes. 
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